Differential Astrometry over 15 



oo 
On 



(N 



ON 



X 



M.A. Perez-Torres a J.M. Marcaide a J.C. Guirado a E. Ros b 

& Depto. de Astronomia, Universitat de Valencia, E-46100 Burjassot, Spain 
h 'Max- Planck- Institut fur Radioastronomie, Auf dem Hiigel 69, Bonn, Germany 



Abstract 



We observed the pair of radio sources 1150+812 and 1803+784 in November 1993 
with a VLBI array, simultaneously recording at 8.4 and 2.3 GHz. We determined 
the angular separation between the two sources with submilliarcsecond accuracy 
by using differential techniques. This result demonstrates the feasibility of high 
■ precision differential astrometry for radio sources separated in the sky by almost 

Q . 15°, and opens the avenue to its application to larger samples of radio sources. 
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1 Introduction 



Very Long Baseline Interferometry (VLBI) is a powerful astrometric technique. 
Centimeter-wavelength VLBI group-delay astrometry of extragalactic radio 
sources routinely provides precisions at the milliarcsecond (mas) level, thus 
allowing a celestial reference frame to be built with milliarcsecond accuracy 
(Ma et al. 1990). One step further in high precision astrometry is phase-delay 
differential astrometry, in which the differences between the phase delays of 
two radio sources are used to determine their relative separation. For close 
source pairs, this technique yields accuracies of a few microarcseconds (//as), 
i.e. the case of the double quasar 1038+528 A and B (Marcaide & Shapiro, 
1983), whose components are separated by less than an arcminute. For sources 
separated up to a few degrees, accuracies of about 0.1-0.3 mas are achieved 
(Guirado et al. 1995, Lara et al. 1996, Ros et al. 1999a). The basic idea un- 
derlying differential astrometry is that the differenced phase-delays from the 
source pair will be largely free of unmodelled systematic effects. If the sources 
are very close, cancellation of systematic effects is almost complete, whereas 
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if the sources lie at increasingly angular distances, the cancellation is less per- 
fect, and the determination of their angular separation is correspondingly less 
accurate. 

In this contribution, we address the question of the applicability of the dif- 
ferential astrometric technique to a pair of strong radio sources, 1150+812 
and 1803+784, separated by almost 15°. We demonstrate the feasibility of the 
technique for this pair of radio sources, and measure their relative separation 
with submilliarcsecond accuracy. In addition, we show that the standard errors 
are dominated by uncertainties in the coordinates of the reference source. 



2 Observations 

We observed the radio sources 1150+812 (QSO) and 1803+784 (BL Lac) on 
November 1993 during 12 hours, simultaneously recording at 8.4 GHz and 
2.3 GHz. We used the following antennas: Effelsberg, Medicina, Onsala, Fort 
Davis, Hancock, Los Alamos, North Liberty, Owens Valley, and the Phased 
VLA. The data were correlated at the Max-Planck-Institut fur Radioastrono- 
mie, Bonn, Germany. We used an observing cycle of 7 min, consisting of 2 
min observing 1150+812, 1.5 min of antenna slew, 2 min observing 1803+784, 
and 1.5 min slew back to 1150+812. Since both sources were strong (> 1 Jy) 
at the epoch of observation, we have detected the radio sources for all of the 
baselines, with high signal-to-noise-ratio (SNR) at both frequencies. 

After correlation, we exported the data and constructed the interferometric 
visibilities using the calibration information provided by the staff at the ob- 
serving antennas. We used the software package DIFMAP (Shepherd et al., 
1995) to obtain images of both sources at each frequency. The images (Figs. 1 
and 2) clearly show the typical core-jet structure for this kind of radio sources, 
the brightest spot being usually associated with the flat-spectrum core, and 
the others with steeper spectrum components of a jet. 1150+812 displays a 
southeastward jet, while 1803+784 displays a westward jet. 



3 Astrometric Data Reduction 

The goal of the data reduction is to obtain a data set of differenced phase- 
delays on which a weighted-least-squares algorithm can be applied to deter- 
mine the relative position of the sources. To obtain this data set, we went 
through several steps, namely, phase connection (or correction of the 27r ambi- 
guities of the phases; Shapiro et al., 1979), structure correction, and removal 
of opacity and propagation medium effects. In the astrometric process, we 
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Fig. 1. Hybrid maps of 1150+812 at 8.4 and 2.3 GHz. Contours are 
0.5,1,2,4,8,16,32,64 and 90% of the peak of brightness, 0.90 Jy/beam and 0.74 
Jy/beam for 8.4 and 2.3 GHz, respectively. The sizes of the natural restoring beams 
are 0.61 x 0.55 mas (PA=15°) at 8.4 GHz, and 2.33 x 1.67 mas (PA=-56°) at 2.3 
GHz. East is left, and North is up. Note that scales are different in both maps. 

used data from the antennas with best performance (Effelsberg, Medicina, 
Fort Davis, Los Alamos, North Liberty, and Owens Valley). 

We used our hybrid maps to correct for the source structure contribution to 
the phase delays, using as reference point the peak of the brightness distri- 
bution of each map. Most of the propagation medium effects are caused by 
the troposphere and the ionosphere. We accurately modelled the tropospheric 
behaviour from detailed information of the weather parameters at each of the 
antenna sites. The ionosphere is a refractive plasma that introduces a dis- 
persive contribution to the phase delays, which scales as v~ 2 . Thus, we used 
our dual-frequency data to estimate and remove such contribution. For it, the 
reference points on the 2.3 GHz maps should correspond to those chosen on 
the 8.4 GHz maps. We took into account opacity effects to obtain a reliable 
map registration (Guirado et al., 1998). 

In this way, we obtained a set of (undifferenced) phase delays, free of struc- 
ture and propagation medium effects. However, the existence of unmodelled 
systematic effects was evident from a visual inspection of the data. We formed 
a set of differenced phase delays by subtracting the delay of each observa- 
tion of 1150+812 from the delay of the previous observation of 1803+784. 
Although the sources are separated by 15° , this procedure largely cancelled 
out systematic effects still present in the undifferenced phase delays. 
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Fig. 2. Hybrid maps of 1803+784 at 8.4 and 2.3 GHz. Contours are 
0.5,1,2,4,8,16,32,64, and 90% of the peak of brightness, 1.48 Jy/beam and 1.35 
Jy/beam for 8.4 and 2.3 GHz, respectively. The size of the natural restoring beam 
is 0.52 x 0.49 mas (PA=61°) at 8.4 GHz, and 2.21 x 1.68 mas (PA=29°) at 2.3 
GHz. East is left, and North is up. Note that scales are different in both maps. 



4 Relative position of the radio sources 



Our estimates of the relative positions of 1150+812 and 1803+784, in right 
ascension and declination, are based on a weighted-least-squares analysis of 
the differenced phase delays. We also included the (undifferenced) phase delays 
for 1803+784 in the analysis to estimate the behaviour of the station clocks. 
The standard deviations of the differenced phase delays, and of the 1803+784 
phase delays were scaled separately so that, for each baseline, the root-mean- 
square of the postfit residuals was unity. The difference between our estimates 
and those given by IERS (1996) are shown in Table 1. 

Table 1 

Estimates of the angular separation of 1150+812 and 1803+784 

Aa - 6 /l 7 m 33f 18473 = -Of 00015 ± Of 00012 
AS + 2°30'25?13601 = 0"00040 ± 0"00030 

The error bars include formal uncertainties from a covariance analysis, and 
residual errors of the geometry of the interferometric array, source position, 
source opacity, and atmosphere. The largest contributor to the error budget in 
the relative separation of the sources is the uncertainty in the a priori position 
of the reference radio source, 1803+784 (IERS, 1996). 
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5 The separation arclength 



Since the angular separation between the two sources is very large, the use 
of the arclength, A, between the sources could be a better estimate of their 
relative separation than Aa and AS, since A is invariant to rotations and, 
therefore, independent of the reference system where the source coordinates 
are defined. The arclength between the two radio sources, compared to the 
arclength as given by IERS is 



A - 14°50'21". 14759 = 0".00040 ± 0".00019 



The uncertainty of the arclength estimate is smaller than for Act and AS, 
as expected. However, the special geometry of the radio sources (see below) 
prevents it from being even smaller. To better clarify this point, we develop 
the expression for the uncertainty in the determination of the arclength, 04, 
using an error propagation equation (Bevington & Robinson, 1992), from 



A(a\, ckb, Sa, Sb) = arccos[cos Sa cos Sb cos(q;b — «a) + sin Sa sin 5 B ] 



where <j\ ~ a Xi a Xj p(xi,Xj) is the covariance of coordinates Xi and Xj (each 
one taking the values a a, cub, Sa, Sb), cr Xi are the a priori uncertainties of the 
source coordinates, and p{xi,Xj) are the correlation coefficients. The North 
Celestial Pole, and the sources 1150+812 (source A) and 1803+784 (source B, 
taken as the reference source), form a triangle of ~ 90° at the North Celestial 
Pole. The geometry is such that cos(a B — «a) <C 1; cos Sa, cos<5 b <sC 1. We can 
then write Eq. 1 approximately as 



(0".aMb) 



bJ ' 



I K 2 cos 2 Sa + cos 2 <5 B 
cos 2 Sa + cos 2 Sb 



where K = <j Sa /o- Sb = 1.4. This gives a A ss 1.18 ag B 0.17 mas. (Note that 
K=l would imply u A rts os B ). This result shows that ~ 89 % of the overall 
standard deviation of the arclength is produced solely by the contribution of 
the a priori standard deviations of the source coordinates. We emphasize that 
this is a consequence of the location of the radio sources in the celestial sphere, 
and dominate by far the standard deviation of our astrometric determination. 
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6 Results and Conclusions 



We observed the pair of strong radio sources 1150+812 and 1803+784 with an 
interferometric VLBI array on November 1993. The antennas simultaneously 
recorded data at 8.4 and 2.3 GHz, which allowed us to remove ionospheric con- 
tributions to the delays. We also estimated the contributions to the delays due 
to the troposphere, source structures, and source opacities. From geometric 
considerations, we infer that the main contribution to the uncertainty in the 
determination of the relative angular separation comes from the uncertainty 
in the a priori position of the reference source. We show that this uncertainty 
is further increased by the particular geometry of the radio source pair. 

We determined with submilliarcsecond accuracy the angular separation of the 
pair of radio sources, separated by almost 15°, and thus demonstrated the 
feasibility of using phase delays from dual-frequency VLBI measurements for 
very large angular separations. This accomplishment opens the applicability 
of the technique to a large sample of radio sources (Ros et al., 1999b). A final 
goal of building a quasi inertial celestial reference frame based on stationary 
radio source cores, with submilliarcsecond accuracy, appears within reach. 



References 

Bevington, P.R. & Robinson, D.K., "Data Reduction and Error Analysis for 

the Physical Sciences", 2nd Ed., McGraw Hill, New York, 1992 
Guirado, J.C., Marcaide, J.M., Elosegui, P., Ratner, M.I., Shapiro, I.I., Eckart, 

A., Quirrenbach, A., Schalinski, C.J., & Witzel, A, 1995, A&A, 293, 613 
Guirado, J.C., Marcaide, Ros, E., Ratner, M.I., Shapiro, I.I., Quirrenbach, & 

Witzel, A, 1998, A&A, 336, 385 
IERS, IERS Annual Report 1995, Observatoire de Paris (1996) 
Lara, L., Marcaide, J.M., Alberdi, A., & Guirado, J.C., 1996, A&A, 314, 672 
Ma, C, Shaffer, D.B., de Vegt. C. Johnston, K.J., & Russell, J.L., 1990, AJ, 

99, 1284 

Marcaide, J.M., & Shapiro, I.I., 1983, AJ, 88, 1133 

Ros, E., Marcaide, J.M., Guirado, J.C., Ratner, M.I., Shapiro, I.I., Krichbaum, 
T.P, Witzel, A., & Preston, R.A., 1999a, A&A, in press. 

Ros, E., Marcaide, J.M., Guirado, & Perez- Torres, M.A., 1999b, A&SS, Pro- 
ceedings of the 3rd Scientific Meeting of the Spanish Astronomical Society. 
Eds. Gorgas, J. & Zamorano, J., in press 

Shapiro, I.I., Wittels, J.J., Counselman III, C.C., Robertson, D.S., Whitney, 
A.R., Hinteregger, H.F., Knight, C.A., Rogers, A.E.E., Clark, T.A., Hutton, 
L.A., & Niell, A.E., 1979, AJ, 84, 1459 

Shepherd, M.C., Pearson, T.J., & Taylor, G.B. 1995, BAAS, 26, 987 



6 



